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Abstract. We have calculated mass- loss rates for a grid 
of wind models covering a wide range of stellar parameters 
and have derived a mass-loss recipe for two ranges of ef- 
fective temperature at either side of the bi-stability jump 
around spectral type Bl. 

For a large sample of O stars, it is shown that there is 
now good agreement between these new theoretical mass- 
loss rates that take multiple scattering into account and 
observations. 

Agreement between the observed and new theoretical 
wind momenta increases confidence in the possibility to 
derive distances to luminous stars in distant stellar sys- 
tems using the Wind momentum Luminosity Relation. 

For the winds of the B stars there is an inconsistency in 
the literature between various mass-loss rate determina- 
tions from observations by different methods. One group 
1 of M determinations of B stars does follow the new theo- 
retical relation, while another group does not. The lack of 
agreement between the observed mass-loss rates derived 
. by different methods may point to systematic errors in 
mass-loss determinations from observations for B stars. 

We show that our theoretical mass-loss recipe is reli- 
able and recommend it be used in evolutionary calcula- 
tions. 

Key words: Stars: early-type - Stars: mass-loss - Stars: 
supergiants - Stars: winds - Stars: evolution 



1. Introduction 

In this paper we present new theoretical mass-loss rates M 
for a wide range of parameters for galactic O and B stars, 
taking multiple scattering into account. These predictions 
for M are compared with observations. The goal of the 
paper is to derive an accurate description of mass loss as 
a function of stellar parameters. 

Early-type stars have high mass-loss rates, which sub- 
stantially affects their evolution (e.g. Chiosi & Maeder 
1986). The winds of early-type stars are most likely driven 
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by radiation pressure in lines and in the continuum. The 
radiation-driven wind theory was first developed by Lucy 
& Solomon (1970) and Castor et al. (1975) (hereafter 
C AK) . At a later stage the theory was improved by Abbott 
(1982), Friend & Abbott (1986), Pauldrach et al. (1986) 
and Kudritzki et al. (1989). 

During the last decade, the radiation-driven wind the- 
ory has been compared with the most reliable mass-loss 
determinations from observations that are available: mass 
loss determined from radio data and from the analysis 
of Ha line profiles. Both Lamers & Leitherer (1993) and 
Puis et al. (1996) came to the conclusion that the theory 
of radiation-driven winds shows a systematic discrepancy 
with the observations. For O stars the radiation-driven 
wind theory predicts systematically lower values for mass 
loss than have been derived from observations. 

Since the discrepancy increases as a function of wind 
density, it is possible that the reason for this is an in- 
adequate treatment of "multiple scattering" in the cur- 
rent state of radiation-driven wind theory. It has been 
suggested (e.g. by Lamers & Leitherer 1993) that the 
"momentum-problem" that has been observed in the 
dense winds of Wolf-Rayet stars is the more extreme ap- 
pearance of this discrepancy seen in the winds of normal 
O-type stars. 

Because the observed mass-loss rates for O type super- 
giants are typically a factor of two higher than the values 
predicted by radiation-driven wind theory, evolutionary 
models would be significantly affected if theoretical values 
were adopted. It is obvious that an accurate description 
of mass loss is of great importance to construct reliable 
evolutionary tracks for massive stars. 

Abbott & Lucy (1985) and Puis (1987) have investi- 
gated the importance of "multiple scattering" relative to 
"single scattering" for the winds of O stars. Abbott & 
Lucy found an increase in M of a factor of about three for 
the wind of the O supergiant £ Puppis if multiple scatter- 
ing was applied in a Monte Carlo simulation. 

We will use a similar Monte Carlo technique in which 
multiple scatterings are taken into account to calculate 
mass-loss rates for a wide range of stellar parameters 
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throughout the upper part of the Hertzsprung-Russell Di- 
agram (HRD). In Sect. |^, the approach to calculate mass- 
loss rates will be briefly described, while in Sect. ||[ a grid 
of wind models and mass-loss rates will be presented. A 
clear separation of the HRD into two parts will be made. 
The first range is that on the "hot" side of the bi-stability 
jump near spectral type Bl, where the ratio of the termi- 
nal velocity to the effective escape velocity at the stellar 
surface (uoo/w CS c) is about 2.6; the second range is that on 
the "cool" side of the jump where the ratio suddenly drops 
to about 1.3 (Lamers et al. 1995). At the jump the mass- 
loss rate is predicted to change dramatically due to a dras- 
tic change in the ionization of the wind (Vink et al. 1999). 
In Sect. U the theoretical wind momentum will be studied 
and in Sect. || fitting formulae for the mass-loss rate will 
be derived by means of multiple linear regression methods: 
this yields a recipe to predict M as a function of stellar 
parameters. In Sect. || these predicted mass-loss rates will 
be compared with observational rates. We will show that 
for O stars theory and observations agree if "multiple scat- 
tering" is properly taken into account. Finally, in Sects. ^| 
and H the study will be discussed and summarized. 

2. Method to calculate M 

The basic physical properties of the adopted Monte Carlo 
(MC) method to predict mass-loss rates are similar to 
the technique introduced by Abbott & Lucy (1985). The 
precise method was extensively described in Vink et al. 
(1999). The core of the approach is that the total loss of 
radiative energy is linked to the total gain of momentum of 
the outflowing material. The momentum deposition in the 
wind is calculated by following the fate of a large number 
of photons that are released from below the photosphere. 

The calculation of mass loss by this method requires 
the input of a model atmosphere, before the radiative ac- 
celeration and M can be calculated. The model atmo- 
spheres used for this study are calculated with the most 
recent version of the non-LTE unified^] Improved Sobolev 
Approximation code (isa-wind) for stars with extended 
atmospheres. For details we refer the reader to de Koter 
et al. (1993,1997). The chemical species that are explic- 
itly calculated in non-LTE are H, He, C, N, O and Si. The 
iron-group elements, which are important for the radiative 
acceleration and M, are treated in the modified nebular 
approximation (see Schmutz 1991). 

3. The predicted mass-loss rates 

Using the procedure summarized in Sect. |[ we have cal- 
culated mass- loss rates for 12 values of T e g in the range 
between 12 500 and 50 000 K. For every effective temper- 
ature a grid of 12 series of models for galactic stars was 
calculated with luminosities in the range log (L*/Lq) = 

1 ISA- wind treats the photosphere and wind in a unified man- 
ner. This is distinct from the so-called "core-halo" approaches. 



Table 1. Parameters for the 12 (L^M*) series of wind 
models. For details about the model assumptions and grid 
spacing, see text. 
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logL. 
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r c 
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(kK) 
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15 


0.055 


14.2 


12.5 - 50.0 


1.3 - 2.6 


2 




20 


0.041 


19.2 


12.5 - 50.0 


1.3 - 2.6 
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5.0 


20 


0.130 


17.4 


12.5 - 50.0 


1.3 - 2.6 


4 




30 


0.087 


27.4 


12.5 - 50.0 


1.3 - 2.6 


5 




40 


0.068 
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12.5 - 50.0 


1.3 - 2.6 
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30 


0.274 


21.8 


12.5 - 50.0 


1.3 - 2.6 


7 




40 


0.206 


31.8 


12.5 - 50.0 


1.3 - 2.6 


8 




50 


0.165 


41.8 
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1.3 - 2.6 
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0.325 


30.4 


12.5 - 50.0 
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Fig. 1. Coverage of the calculated wind models over the 
HRD. The crosses indicate the model values of log L/Lq 
and T ff- Evolutionary tracks from Meynet et al. (1994) 
are shown for several initial masses, which are indicated 
in the plot. The Zero Age Main Sequence (ZAMS) is also 
shown. 



4.5 - 6.25 and masses in the range = 15 - 120 M@. 
For these 144 models, mass-loss rates were calculated for 
three values of the ratio Woo/w CS c, yielding a total number 
of 432 models. 

The parameters for all series of models are indicated 
in Table |l|. In Fig. [j] the luminosities and effective temper- 
atures of the models are indicated with asterisks on top of 
evolutionary tracks to show the coverage of the model grid 
over the upper HRD. To study the mass-loss dependence 
on different stellar parameters (L, M and T e s) separately, 
a wide range of parameters was chosen, this implies that 
some of the models in Fig. [I] have positions to the left of 
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the main sequence. We enumerate the assumptions in the 
model grid: 

1. The models are calculated for solar metallicities (Allen 
1973). 

2. The stellar masses in the grid of models were chosen 
in such a way that they are representative for the evo- 
lutionary luminosities of the tracks from the Geneva 
group (Meynet et al. 1994). To investigate the depen- 
dence of M on Mi,, a number of smaller and larger 
values for M* was also chosen (see column (3) in Ta- 
ble@). 

3. The grid was constructed in a way that T e ^ 0.5 (see 
column (4) in Table ^) , where T e is the ratio between 
the gravitational acceleration and the radiative accel- 
eration due to electron scattering. T e is given by: 

La P - f T. \ / M \ - 1 

J- e 



Table 2. Bi-stability jump characteristics for the 12 
(L*, M») series of wind models. 



= 7.66 10- J cr e (1) 

AwcGM e \L G J \MqJ 1 ' 

where <j e is the electron scattering cross-section (its 
value is taken as determined in Lamers & Leitherer 
1993) and the other constants have their usual mean- 
ing. For values of T e > 0.5, the stars approach their 
Eddington limit and the winds show more dramatic 
mass-loss behaviour. In this study, stellar parameters 
for these "Luminous Blue Variable-like" stars are ex- 
cluded to avoid confusion between various physical 
wind effects. 

All series of models from Table |l| have effective tem- 
peratures between 12 500 and 50 000 K, with a stepsize 
of 2 500 K from 12 500 to 30 000 K and a stepsize of 
5 000 K, starting from 30 000 up to 50 000 K. 
We calculated M for wind models with a /3-type veloc- 
ity law for the accelerating part of the wind: 



v(r) 



1 



R* 







(2) 



Below the sonic point, a smooth transition from this 
velocity structure is made to a the velocity that follows 
from the photospheric density structure. A value of 
/3 = 1 was adopted in the accelerating part of the wind. 
This is a typical value for normal OB supergiants (see 
Groenewegen & Lamers 1989; Haser et al. 1995; Puis 
et al. 1996). At a later stage models for other (3 values 
will be calculated and it will be demonstrated that the 
predicted M is essentially insensitive to the adopted 



value of f3 (see Sect. p.4[ ). 

The dependence of M on various values of i>oo was de- 
termined. Lamers et al. (1995) found that the ratio 
Vao/vcsc ~ 2.6 for stars of types earlier than Bl, and 
drops to Voo/v csc ~ 1.3 for stars later than type Bl. 
Therefore, we calculated mass-loss rates for various in- 
put values of this ratio, namely 1.3, 2.0 and 2.6 to in- 
vestigate the mass loss as a function of this parameter, 
similar to that in Vink et al. (1999). For the determi- 
nation of v csc , the effective mass M e g = M*(l — T e ) 
was used. M e s is given in column (5) of Table @. 



series 


log L, 


M* 


log(AM) 


T-ijump 


< p > jump 


no 


(£©) 


(Mo) 




(K) 


(g cm -3 ) 


1 


4.5 


15 


0.78 


23 750 


-14.82 


2 




20 


0.61 


22 500 


-15.13 


3 


5.0 


20 


0.83 


26 250 


-14.22 


4 




30 


0.87 


25 000 


-14.68 


5 




40 


0.73 


25 000 


-14.74 


6 


5.5 


30 


0.76 


26 250 


-13.89 


7 




40 


0.81 


26 250 


-14.13 


8 




50 


0.82 


25 000 


-14.40 


9 


5.75 


45 


0.77 


25 000 


-13.93 


10 


6.0 


60 


0.76 


25 000 


-13.66 


11 




80 


0.76 


26 250 


-13.89 


12 


6.25 


120 


0.77 


25 000 


-13.87 



3.1. M for supergiants in Range 1 
(30 000 < T efi < 50 000 K) 

The results for the complete grid of all the 12 (L„, M*) se- 
ries are plotted in the individual panels of Fig. ||. Note that 
for each calculated point in the grid, several wind models 
had to be calculated to check which adopted mass-loss rate 
was consistent with the radiative acceleration (see Lucy & 
Abbott 1993). This yields predicted, self-consistent values 
for M (see Vink et al. 1999). 

For each (L*, M*) set and for each value of Voo/v csc , we 
found that the mass loss decreases for decreasing effective 
temperature between 50 000 and 27 500 K. The reason 
for this fall-off is essentially that the maximum of the flux 
distribution gradually shifts to longer wavelengths. Since 
there are significantly less lines at roughly A <; 1800 A 
than at shorter wavelength, the line acceleration becomes 
less effective at lower T c g, and thus M decreases. 

3.2. M at the bi-stability jump around 25 000 K 

Between about T e g = 27 500 and 22 500 K the situation is 
reversed: in this range the mass loss increases drastically 
with decreasing T c g. These increments in M coincide both 
in T c ff and in size of the M jump with the bi-stability jump 
that was studied by Vink et al. (1999). They showed that 
the origin of the jump is linked to a shift in the ionization 
balance of iron in the lower part of the wind and that it is 
this element that dominates the line driving at the base of 
the wind. Below T cS ~ 25 000 K, Fe iv recombines to Fe 
in and as this latter ion is a more efficient line driver than 
Fe iv, the line acceleration in the lower part of the wind 
increases. This results in an upward jump in M of about a 
factor of five and subsequently a drop in v^. The drop in 
Vqo was predicted to be a factor of two, which is confirmed 
by determinations of of OB supergiants from ultravi- 
olet data by Lamers et al. (1995). A comparison between 
the spectral type of the observed bi-stability jump and the 
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Fig. 2. The calculated mass-loss rates M as a function of T c s for the grid of 12 (L*, M„) series for three values of the 
ratio Woo/vcsc- The values for Voo/?W are indicated in the lower part of the first panel. The stellar parameters L* and 
M* are indicated in the upper part of each panel. The thin dotted lines connect the calculated mass-loss rates. The 
thick solid lines indicate two multiple linear regression fits to the calculated values (see Sect. |]). Note that some of the 
panels seem to indicate that a more accurate fit may be possible. However, the fits have been derived by multiple linear 
regression methods and thus the thick solid lines show the mass loss as a function of more than just one parameter. 
All models were calculated for solar metallicities. 
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Fig. 3. Characteristic < p > and T e g of the bi-stability 
jump around T e g = 25 000 K. An explanation for the 
different symbols is given in the legend of the plot. The 
solid line represents a linear fit through the average jump 
parameters log < p > and T e g. 

effective temperature of the predicted jump, was made in 
Vink et al. (1999). 

Since we know from both theory and observations that 
the ratio Voo/v esc jumps from ~ 2.6 at the hot side of 
the jump to ~ 1.3 at the cool side of the jump, we can 
predict the jump in mass loss for all 12 (Z*,M*) series 
of models. The size of the jump is defined as the differ- 
ence between the minimum M at the hot side of the jump 
(where foo/^csc = 2.6) and the maximum M at the cool 
side (where Voo/w sc — 1-3) in Fig. ^. The size of the pre- 
dicted jump in M (log AM) is indicated in column (4) of 
Table |^: AM is about a factor of five to seven. Table | 
tabulates additional characteristics for the models at the 
bi-stability jump. 

The jump in mass loss around T e g — 25 000 K is not 
exactly the same for all series of models: the jump occurs 
at somewhat different effective temperatures. This is no 
surprise, since the ionization equilibrium does not only 
depend on temperature, but on density as well. A smaller 
value of the ratio Woo/w eS c leads to a larger density in the 
wind. Hence, the jump is expected to start at higher T e g 
for smaller Voo/fosc- This behaviour for the position of T e g 
of the jump is confirmed by all individual panels in Fig. 0. 
To understand the behaviour of the bi-stability jump as a 
function of the other stellar parameters, i.e. M* and L*, 
we will compare the wind characteristics of the 12 series of 
models around the bi-stability jump in some more detail. 

First we define a characteristic wind density at 50 % 
of the terminal velocity Vqo of the wind: < p >. For a 
standard velocity law with (3 = 1, this characteristic wind 
density is given by 

M 

<P> = a m (3) 



For all 12 series of models this characteristic density < p > 
is plotted vs. the effective temperature of the jump. This 
is done for both the minimum M (at the hot side of the 
jump) and the maximum M (at the cool side of the jump). 
Figure [| shows the location of the bi-stability jump in 
terms of T e g as a function of < p >. The characteristic 
densities and effective temperatures for the cool side of 
the jump are indicated with "diamond" signs and with 
"plus" signs for the hot side. As expected, for all 12 mod- 
els the minimum M corresponds to a relatively low p and 
relatively high T c g, whereas the maximum M corresponds 
to a relatively high p, but low T c g. Note that the effective 
temperature at minimum and maximum mass loss is not 
a very smooth function of wind density. This is due to 
our choice of resolution in effective temperature of the 
grid. We have checked whether the obtained minima and 
maxima were indeed the extreme mass-loss values by cal- 
culating extra models at intermediate values of T e g. The 
minimum and maximum M values obtained with the ini- 
tial grid resolution were found to be similar to those de- 
termined with a the finer resolution. We thus concluded 
that the initial resolution of the grid was justified. 

The "filled circles" represent the average values of T c g 
and < p > for the "jump" model for each M*) series. 
The "jump" model is a hypothetical model between the 
two models where M is maximal and minimal. The solid 
line indicates the best linear fit through these averages. 
The relation between the jump temperature (in kK) and 
log < p > is given by: 

T jum P = 49>1 £ ± g ^ + Lg7 £ ± 64 j i og < p > (4) 

The average temperature and density of the jump are 
given in columns (5) and (6) of Table ||. Note that the 
range in T^ mp is relatively small; all 12 series of mod- 
els have jump temperatures in the range between 22.5 
£ T cS <, 26 kK. 

Figure [| shows the behaviour of the characteristic den- 
sity log < p > as a function of T e . Again this is done for 
both the cool and hot side of the jump, and for the av- 
erage between them. As expected, log < p > increases as 
F e increases. Since the average characteristic wind density 
at the jump shows an almost linear dependence on r e , a 
linear fit through the average densities is plotted. This is 
the solid line in Fig. [|. The relation between log < p > 
and r e is given by: 

log < p > = - 14.94 (± 0.54) + 3.2 (± 2.2) T e (5) 

From the quantities L* and M* it is now possible to esti- 
mate log < p > using Eq. (||) and subsequently to predict 
X^ mp using Eq. (Q). Later on this will be used as a tool 
to connect two fitting formulae for the two ranges in T e g 
at either side of the bi-stability jump (see Sect. ||). 

3.3. M for supergiants in Range 2 
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ture the mass-loss rate still decreases and for which mod- 
els approaching the second bi-stability jump, the mass loss 
again increases. 

The relation found between the temperature of the sec- 
ond bi-stability jump and log < p > is determined by eye 
and is roughly given by: 



T jump2 = ioo + 6 log < p > 



(6) 



where T> ump2 is in kK. From the quantities and M* 
it is again possible to estimate log < p > using Eq. (JsJ) 
and then to roughly predict T> ump2 using Eq. (||). This 
formula will be used for our mass-loss recipe at the low 
temperature side (see Sect. B). 



Fig. 4. Characteristic < p > at the bi-stability jump as a 
function of T e . An explanation for the different symbols 
is given in the legend of the plot. The solid line indicates 
a linear fit through the average jump parameters for log 
< p >. 



(12 500 < T cff < 22 500 K) 

Figure H shows that at effective temperatures T e g < 22 500 
K, M initially decreases. This is similar to the M be- 
haviour in the T e g range between 50 000 and 27 500 K. 
For some series (dependent on the adopted L^/M*) the 
mass loss decreases until our calculations end at T e g = 
12 500. For other series of L* and Af*, the initial decrease 
suddenly switches to another increase. Vink et al. (1999) 
already anticipated that somewhere, at lower T e fj, a re- 
combination would occur from Fe in to n similar to the 
recombination from Fe iv to ill at ~ 25 000 K. Lamers et 
al. (1995) already mentioned the possible existence of such 
a second bi-stability jump around T c g — 10 000 K from 
their determinations of Voc/v CSCl but the observational ev- 
idence for this second jump is still quite meagre. 

3.4- M at the second bi-stability jump around 12 500 K 

To understand the characteristics of the "second" bi- 
stability jump as a function of different stellar parameters 
(M* and L*), we have also studied the models around this 
second jump in some more detail. 

Since our model grid is terminated at 12 500 K, it is 
not possible to determine the maximum M of the second 
bi-stability jump in a consistent way, similar to that of 



the first jump discussed in Sect. 3.2. Thus, it is not possi 



ble to determine the exact size of the second jump in M . 
Neither is it possible to derive an accurate equation for 
the position of the second bi-stability jump in T e g (as was 
done in Eq. || for the first jump around 25 000 K). Still, 
it is useful to determine a rough relationship between the 
position of the second jump in T e g and the average log 
< p > by investigating for each model at which tempera- 



4. The wind momentum 

4--1- The wind efficiency number n 

In this section, we present values for the wind efficiency 
number rj for the different (L*,M*) series, r\ (sometimes 
called the wind performance number) describes the frac- 
tion of the momentum of the radiation that is transferred 
to the ions in the wind: 



(7) 



Figure |5| shows the behaviour of ry as a function of T e g for 
the complete grid of models. Figure |5| demonstrates that 
rj is not constant as a function of T c g. The figure shows 
that when a star evolves redwards at constant luminosity 
(from high to low temperature) the momentum efficiency r\ 
initially decreases until the star approaches the bi-stability 
jump around 25 000 K, where the wind efficiency suddenly 
increases by a factor of two to three. Subsequently, below 
about 22 500 K, rj decreases again and in some cases (again 
dependent on L* and Af*) it eventually jumps again at 
the second bi-stability jump. This overall behaviour of rj 
is similar to that of M as shown in Fig. ^. 

In some of the panels of Fig. 0, i.e. in those cases where 
L*/M* is large, rj exceeds the single scattering limit. 



n 



L*/c 



> 1 



(8) 



This occurs at T cS > 40 000 K and log {L*/L & ) > 6. It 
suggests that already for high luminosity OB stars stellar 
winds cannot be treated in the single scattering formalism. 
The single-scattering limit which is definitely invalid for 
the optically thick winds of Wolf-Rayet type stars, is often 
assumed to be valid for the winds of "normal" supergiants. 
Here, however, we come to the conclusion that due to mul- 
tiple scattering, n already exceeds unity for luminous, but 
"normal" OB supergiants, in case log(L/L Q ) <; 6. This 
was already suggested by Lamers & Leitherer (1993) on 
the basis of observations. 
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4-2. The importance of multiple scattering 

Puis et al. (1996) proposed that the reason for the sys- 
tematic discrepancy between the observed mass-loss rates 
and recent standard radiation driven wind models (Paul- 
drach et al. 1994) was caused by an inadequate treatment 
of multi-line effects in these wind models. To compare 
our new mass-loss predictions with the most sophisticated 
prior investigations, it is useful to briefly discuss the most 
important assumptions that are made in modelling the 
wind dynamics of OB- type stars. The following four basic 
choices must be made: 

1. One may treat the photosphere and wind in a "core- 
halo" approximation, or one may not make this dis- 
tinction and treat photosphere and wind in a "unified" 
way. This choice must be made twice, i.e. with respect 
to the calculation of the occupation numbers as well 
as with respect to the computation of the line force. 

2. One may adopt a "single-line" approach, i.e. neglecting 
effects caused by overlapping lines, or one may follow 
an approach including "multi-line" effects. 

3. One solves the rate equations for all relevant ions ex- 
plicitly in non-LTE, or one adopts a "nebular type of 
approach" to calculate the ionization balance. 

4. One solves the equation of motion self-consistently, or 
one derives the wind properties from a global energy 
argument. 

Standard radiation driven wind models (CAK, Abbott 
1982, Pauldrach et al. 1994) treat the momentum equation 
in a core-halo approach (1) adopting the single- line ap- 
proximation (2). Various degrees of sophistication can be 
applied to determine the occupation numbers. The studies 
of Pauldrach et al. (1994) and Taresch et al. (1997) repre- 
sent the current state-of-the-art, i.e. they treat all relevant 
ions explicitly in non-LTE (3) and solve the equation of 
motion self-consistently (4). Pauldrach et al. (1994) also 
use a unified method for the calculation of the occupa- 
tion numbers, but a "core-halo" approach is applied with 
respect to the line force. Additionally, as line overlap is 
neglected in the method used by Pauldrach et al. (1994), 
these models can overestimate the line force as unattenu- 
ated photospheric flux is offered to each line, which conse- 
quently may produce efficiency numbers larger than unity. 

Puis (1987) found that for winds of relatively low den- 
sity (say n 1/2) the inclusion of multi-line effects leads 
to a reduction of wind momentum compared to the stan- 
dard model due to backscattering and blocking of photons 
in the lower part of the wind. For winds of relatively high 
density (say n ^ 1), such as the dense winds of Wolf-Rayet 
stars, the situation is likely to be reversed. Here momen- 
tum transfer from an extended diffuse field is expected 
to dominate over the effect of the attenuation of flux in 
the layers just above the photosphere. This could result 
in more mass loss compared to the standard radiation 
driven wind theory (Abbott & Lucy 1985, Springmann 



1994). Wolf-Rayet and Of/WN stars profit from a layered 
ionization structure, which increases the number of lines 
that can be used for the driving and thus increasing the 
mass loss (Lucy & Abbott 1993, de Koter et al. 1997). 

Our method differs in almost all aspects from that of 
Pauldrach et al. (1994). In our method, photosphere and 
wind are treated in a unified manner (1) and we properly 
take multi-scatterings into account with a Monte Carlo 
technique (2). On the other hand, we derive the level pop- 
ulations of the iron-group elements using (a sophisticated 
version of) the nebular approximation (3). Finally, we de- 
rive the mass loss from a global energy argument (4) . This 
distinct difference of approach implies that a comparison 
between both methods is difficult. Still, we will address 
some of the differences in approach by focusing on a star 
with parameters representative for the 04I(f)-star £ Pup- 
pis, which has been studied in detail by Abbott & Lucy 
(1985), Puis (1987), Pauldrach et al. (1994) and Puis ct 
al. (1996). 

We can test the difference between single scattering 
and multiple scattering by allowing photons to interact 
with a line only once. Fig. ^ show a comparison between 
the single- and multiple scattering case for three repre- 
sentative wind models at T c g = 40,000 K. The model pa- 
rameters arc given in Table 3. For the often studied wind 
of the O supergiant £ Puppis, which has a mass-loss rate 
of M ohs = 5.9 x lO^Afoyr- 1 (Puis et al. 1996), the ob- 
served efficiency number is about r\ ~ 0.6, suggesting that 
the real efficiency of multiple vs. single scattering is a fac- 
tor of about four for £ Puppis (see Fig. ^). This is close 
to the findings of Abbott & Lucy (1985) who found an 
increase in M by a factor of 3.3 for the wind of £ Puppis 
if multiple scattering was taken into account in a Monte 
Carlo simulation. 

Note from the figure that at low wind densities, the 
single- and multiple scattering approach converge, as one 
would expect. For typical O-stars, which have n ^ 0.5, 
the mass loss will increase by up to a factor of two when 
multiple scattering is properly included. The Wolf-Rayet 
stars, located at the extreme high wind density side, and 
which in some cases have observed efficiency numbers of 
factors 10 or even higher, may benefit by factors of up to 
- 50. 

The reason why Puis (1987) found a reduced mass 
loss for £ Puppis when comparing the single-line approach 
with the multi-line approach is because the single-line ap- 
proach (which is not the same as the single scattering pro- 
cess) overestimates the line force at the base of the wind, 
where the mass loss is fixed. However, a similar relative 
behaviour is not found when we compare the predicted 
single-line mass loss M = 5.1 x 10 _6 MQyr _1 of Pauldrach 
et al. (1994) with the value of M = 8.6 x 10" 6 M o yr- 1 
derived from our fitting formula based on multiple scat- 
tering models. It is not possible to exactly pinpoint the 
cause of this difference, but it is likely to be related to 
differences in our multi-line treatment and that of Puis 
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Table 3. The relative importance of multiple (MS) vs. In the theory of line driven winds, the reciprocal value 

single scattering (SS) for a wind model at T c g = 40 000 of x equals (Puis et al. 1996): 

K. 

-1/.t = a e g = a — 5 (10) 



To 


log L* 

(L S ) 


M* 
(M ) 


77 MS 


log M ss 


log M MS 


M MS 

mss 


0.041 


4.5 


20 


0.107 


-7.87 


-7.72 


1.41 


0.206 


5.5 


40 


0.460 


-6.46 


-5.95 


3.24 


0.434 


6.0 


60 


1.07 


-5.76 


-4.97 


6.17 



Here a and 8 are force multiplier parameters, describ- 
ing the radiative line acceleration gn nc through the stellar 
wind: 



Sline OC 



T3 2 



where n c is the electron density and W is the dilution 
factor, a corresponds to the power law exponent of the 
line strength distribution function controlling the relative 
number of strong to weak lines. If only strong (weak) lines 
contribute to the line acceleration force, then a = 1 (0). 
The predicted value of a is about 0.6. The parameter S de- 
scribes the ionization balance of the wind. Values for this 
parameter are usually between 0.0 and 0.1. For a detailed 
discussion of the paramctcrisation of the line acceleration, 
see e.g. CAK, Abbott (1982) and Kudritzki et al. (1989). 

The important point to note here is that possible 
changes in the slope function of effective temper- 

ature reflect the fact that the stellar winds are driven by 
different sets of ions, i.e. lines of different ions. Figure || 
shows that around the bi-stability jump at T e g ~ 25 000 K, 

, . n . r V increases for decreasing T c ff . This implies that one does 

Fig. 6. 1 he efficiency of multiple-scattering for a range of , . , , 

. , , . . . , . , . , „„ not necessarily expect a universal WLR over the complete 

wind densities. Mb refers to multiple-scattering, and bb , c „ „ , . , 

, . . ° spectral range of (J, B and A stars, nor does one expect a 

rclcrs to single-scattering. , £ r j-xr i ^ 

° ° constant value of a c g or x for different spectral types. 



I i i i i i i i i i i i i i i i i i i i i , , 
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7] 



(1987). Contrary to Puis (and also contrary to Abbott & 
Lucy 1985), we do not adopt the core-halo approximation. 
The formation region of the strong driving lines extends 
from the photosphere out to the base of the wind. If one 
assumes an input photospheric spectrum representative 
of the emergent ultraviolet spectrum as in a core-halo ap- 
proach, one may overestimate the blocking in the subsonic 
wind regime which results in a lower mass loss. 

4-3. The Modified Wind Momentum II 

Kudritzki et al. (1995) have defined the Wind momentum 
Luminosity Relation (WLR): 

LI = M Voo R°/° oc LI (9) 

where LI (or MvoqR® 5 ) is called the "modified wind mo- 
mentum" . Observations of M and v^o of O supergiants 
have shown that log LI is proportional to log L* (see e.g. 
Puis et al. 1996). The WLR may in principle be used as 
a tool to derive distances to galaxies (see Kudritzki et al. 
1995). 



5. Mass loss recipe 

In this section we present a theoretical mass loss formula 
for OB stars over the full range in T e g between 50 000 
and 12 500 K. The mass-loss rate as a function of four 
basic parameters will be provided. These parameters are 
the stellar mass and luminosity, effective temperature and 
terminal velocity of the wind. To obtain a mass-loss recipe, 
we have derived interpolation formulae from the grid of 
M calculations presented in Sect. [|. The fitting procedure 
was performed using multiple linear regression methods 
to derive dependence coefficients. We have applied this 
method for the two ranges in T e g separately. The first 
range is roughly the range for the O-type stars between 
T cS = 50 000 and 30 000 K. The second range is between 
T e g = 22 500 and 15 000 K, which is roughly the range for 
the B-type supergiants. The two relations are connected 
at the bi-stability jump. We have already derived the jump 
parameters for different series of models in Sect. |^, so we 
have knowledge about the position of the jump in T e s 
as a function of stellar parameters. This will be applied in 
the determination of mass loss for stars with temperatures 
around the bi-stability jump. 
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5.1. Range 1 (30 000 < T cS < 50 000 K) 

The first range (roughly the range of the O-type stars) is 
taken from T off between 50 000 K and 30 000 K. In this 
range the step size in effective temperature equals 5 000 
K. So, for the first range we have five grid points in T e ff. 
Five times 12 series of (L», M*), together with three ratios 
of ("Uoo/vcsc) yields a total of 180 points in M for the first 
range. We have found that for the dependence of M on 
T e ff, the fit improved if a second order term (log T c s) 2 was 
taken into account. In fact, this is obvious from the shapes 
of the plots in the panels of Fig. |[ The best fit that was 
found by multiple linear regression is: 



1.601 (±0.055) log 



2.0 



log M 



- 6.697 (±0.061) 

+ 2.194 (±0.021) log(L„/10 5 ) 

- 1.313 (±0.046) log(M*/30) 



1.226 (±0.037) log 



2.0 



+ 0.933 (±0.064) log(T off /40000) 
- 10.92 (±0.90) {log(T off /40000)} 2 

for 27 500 < T cff < 50 000K 



(12) 



where M is in M Q yr -1 , and M* are in solar units 
and T ff is in Kelvin. Note that M» is the stellar mass not 
corrected for electron scattering. In this range ^oo/^csc = 
2.6. Equation |l2| predicts the calculated mass- loss rates of 
the 180 models with a root-mean-square (rms) accuracy 
of 0.061 dex. The fits for the various (L*,M*) series are 
indicated with the thick lines in the panels of Fig. Note 
that some of the panels in Fig. || seem to indicate that 
a more accurate fit might have been possible. However, 
Eq. (12) is derived by multiple linear regression methods 
and thus it provides the mass loss as a function of more 
than just one parameter. 

5.2. Range 2 (15 000 < T cff < 22 500 K) 

The second range (roughly the range of the B-type super- 
giants) is taken from T e ff between 22 500 and 15 000 K. 
In this range the step size in effective temperature equals 
2 500 K. For this range, there are four grid points in T e g. 
Four times 12 series of (L*,M»), together with three ra- 
tios of (voo/fcsc) yields a total of 144 points in M. In this 
range the fit did not improve if a second order term in 
effective temperature was taken into account, so this was 
not done. The best fit that was found by multiple linear 
regression for the second range is: 



log M = - 6.688 (±0.080) 

+ 2.210 (±0.031) log(L*/10 5 ) 
- 1.339 (±0.068) log(M*/30) 



1.07 (±0.10) log(T cff /20000) 



for 12 500 < r efT < 22 500K 



(13) 



where again M is in Mq yr -1 , L* and M* are in solar 
units and T e ff is in Kelvin. In this range Voo/fcsc = 1-3. 
The fitting formula is also indicated by solid lines in the 
panels of Fig. ||. Equation [l3] predicts the calculated mass- 
loss rates of the 144 models for this T c g range with an 
rms accuracy of 0.080 dex. For this second range (12 500 
< T cff < 22 500 K) the fit is slightly less good than for the 
first T c ff range. This is due to the presence of the second 
bi-stability jump which already appears in some (L»/M») 
cases, as was shown in Fig. |[ If those models that do show 
the second bi-stability jump, i.e. stars with high r e , are 
omitted from the sample, the accuracy improves to ~ 0.06 
dex. In all cases the rms is Ss 0.08 dex in log M, which im- 
plies that the fitting formulae yield good representations 
of the actual model calculations. 

We are aware of the fact that there could be systematic 
errors in our approach, since we have made assumptions 
in our modelling. For a discussion of these assumptions, 
see Vink et al. (1999). Whether there are still systematic 
errors between the observed mass-loss rates and these new 
predictions of radiation-driven wind theory, will be inves- 
tigated in Sect. |[ 

5.3. The complete mass-loss recipe 

For stars with effective temperatures higher than 27 500 
K, one should apply the mass-loss formula for the first 
range (Eq. |lj) ; for stars with T e g lower than 22 500 K the 
formula for the second range (Eq. |l^) is to be used. In the 
range between 22 500 and 27 500 K, it is not a priori known 
which formula to apply. This due to the presence of the 
bi-stability jump. Nevertheless, it is possible to retrieve a 
reliable mass-loss prediction by using Eqs. [| and || as a 
tool to determine the position of the jump in T e g. 

In predicting the mass-loss rate of stars close to the bi- 
stability jump, one should preferentially use the observed 
^oo/vcsc value to determine the position with respect to 
the jump. This is a better approach than to use the tools 
from Eqs. [| and |^ to determine the position of the jump. 
The reason is that errors in the basic stellar parameters 
may accidcntly place the star at the wrong side of the 
jump. A computer routine to calculate mass loss as a func- 
tion of input parameters is available either upon request 
or at the following url: www. astro. uu.nl/^jvink/. 

If Woo is not available, as is the case for evolutionary 
calculations, one should adopt the ratio Voo/v csc — 2.6 for 
the hot side of the jump and i>oo/v C sc = 1-3 for the cool 
side of the jump, in agreement with the analysis by Lamers 
et al. (1995). Note that the exact T e g of the jump is not 
expected to have a significant effect on evolutionary tracks 
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calculated with this new mass-loss description, since the 
most luminous stars spend only a relatively short time 
around T c g ~ 25 000 K during their evolution. 

Since our calculations were terminated at 12 500 K, 
we are not able to determine the size and the position of 
the second bi-stability jump. Predicting the mass-loss be- 
haviour below this second jump would therefore be spec- 
ulative. Yet, for evolutionary tracks the mass loss below 
12 500 K is an important ingredient in the evolutionary 
calculations. We roughly estimate from our grid calcula- 
tions that for a constant ratio of Woo/fesc the increase in 
M around 12 500 is about a factor of two, similar to that 
found for the first jump near 25 000 K. Furthermore, ob- 
servations by Lamers et al. (1995) indicate that for stars 
around 10 000 K, Voo/vcsc drops again by a factor of two 
from Voo/fcsc — 1-3 to about 0.7. It is therefore plausible 
to expect that the size in M of the second jump is compa- 
rable to the size of the first jump. So, AM of the second 
jump is expected to be a factor of five also. We argue that 
this second jump should also be considered in evolution- 
ary calculations and suggest Eq. (|l^) could be used for 
effective temperatures below the second jump when the 
constant in Eq. (|l3| ) is increased by a factor of five (or log 
AM = 0.70) to a value of -5.99. The mass-loss recipe can 
be applied for evolutionary calculations until the point in 
the HRD where line driven winds become inefficient and 
where probably another mass-loss mechanism switches on 
for the cooler supergiants (see Achmad et al. 1997). We 
suggest that in the temperature range below the second 
jump «oo /f osc = 0.7 is adopted. 

5.4- The dependence of M on the steepness of the 
velocity law [3 

To test the sensitivity of our predictions of mass-loss rates 
on different shapes of the velocity law, we have calculated 
series of models for (3 — 0.7, 1.0 and 1.5. This is a rea- 
sonable range for OB stars, see Groenewegen & Lamers 
1989; Puis et al. 1996). The adopted stellar parameters 
for this test are L* = 10 5 L and M, = 20 M Q . We have 
calculated M for all of the above (3 values for wind models 
with the three values / v osc = 2.6, 2.0 and 1.3. 

From the results shown in Fig. [7] we derived for the 
dependence of M on the adopted value of (3: 

log M = C + 0.112 (±0.048) log(/?/1.0) (14) 

where C is a constant. This relation is valid for the range 
between (3 — 0.7 - 1.5. Since the dependence on this pa- 
rameter is significantly smaller than that on the other 
parameters, L*, M*, T c g and Voo/f OS c, as was found in 
Eqs. [l2] and [l3|, we have omitted the (3 dependence from 
the mass loss recipe. We have just presented the (3 depen- 
dence in this section for the sake of completeness, but we 
can conclude that the predicted mass-loss rates are only 
marginally sensitive to the shape of the adopted velocity 
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Fig. 7. Dependence of M on the shape of the velocity law, 
for three values of (3 = 0.7, 1.0 and 1.5, as is indicated in 
the lower left corner of the plot. The values for Voo/v csc 
are indicated in the upper left corner of the plot. For other 
stellar parameters, see text. 

law. One could argue that a (3 dependence on M could 
be of significance for more extreme series of models. This 
was tested, but it turned out that for the high T e series, 
the (3 dependence is also insignificant, i.e. deviations of 
predicted M are less than AM <S 0.03 dex. This shows 
that we can safely omit the (3 dependence on M in the 
mass-loss recipe for the O and B stars. 

6. Comparison between theoretical and 
observational M 

6.1. M comparison for Range 1 
(27 500 < T off < 50 000 K) 

An extended compilation of observed mass-loss rates of 
early-type OBA stars is obtained by Lamers et al. (2000; 
in preparation). Since both the ultraviolet and the infrared 
method do not yet yield reliable rates, only mass-loss rates 
based on radio free-free emission and emission of Ha have 
been considered. The Ha mass-loss rates and their stellar 
parameters arc from: Herrero et al. (2000); Kudritzki et 
al. (1999); Lamers & Leitherer (1993) (these Ha equiva- 
lent width values are corrected with the curve of growth 
method from Puis et al. 1996); Puis et al. (1996); Scud- 
cri et al. (1992), Scuderi (1994), and Scuderi & Panagia 
(2000). The radio rates are from the compilation of Lamers 
& Leitherer (1993); from Leitherer et al. (1995) and Scud- 
eri et al. (1998). The observed terminal velocities are from 
the same papers. These were mainly determined from P 
Cygni profiles. The stellar masses are derived from evo- 
lutionary tracks of Meynet et al. (1994). For a critical 
discussion of the observed mass-loss rates and for the se- 
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lection of the most reliable data, see Lamers et al. (in 
preparation) . 

For all these stars with known observational mass-loss 
rates and stellar parameters, we have determined theoret- 
ical M values with the mass-loss recipe that was derived 
in Sect. |[ A star-to-star-comparison between these pre- 
dicted mass-loss rates and those derived from observations 
is presented in Fig. |L In this plot only the stars above the 
bi-stability jump (where T c g > 27 500 K) are included. 
The mass-loss rates from Puis et al. (1996) are indicated 
with a different symbol (filled circle), because these are 
obtained from a homogeneous set, and are analyzed with 
the most sophisticated wind models. Note that the out- 
lier at log M bs ^ - 7.4 is the star C Oph (HD 149757) 
for which Lamers & Leitherer (1993) reported that the 
mass-loss rate is uncertain. 

The errors in Fig. || can be due to several effects. There 
is an error in the theoretical fitting formula, though this 
error is only 0.061 dex (see Sect. 5.1). There could also 



be systematic errors due to assumptions in the modelling. 
Furthermore, there could be systematic errors in the mass- 
loss determinations from observations. Such systematic ef- 
fects may for instance occur if the clumping factor in the 
wind changes with distance to the central star. This be- 



cause the Ha and radio emission originate from distinctly 
different regions in the stellar wind. However, Lamers & 
Leitherer (1993) have shown that for a significant sample 
of O stars there is good agreement between the radio and 
the Ha mass-loss rates. 

The random errors in the observational mass-loss rates 
are due to uncertainties in the stellar parameters and in 
the mass-loss determinations. We tentatively estimate the 
intrinsic errors in the observed mass-loss rates from the 
radio and Ha method to be on the order of 0.2 - 0.3 dex 
(see Lamers et al. in preparation). This means that for 
a star-by-star comparison between observations and the- 
ory one would expect a scatter around the mean which 
is a combination of the theoretical and observed uncer- 
tainties. This error is on the order of 0.3 dex. The scatter 
between observational and theoretical mass-loss rates for 
the O stars from Fig. ^ that was actually derived, equals 
0.33 dex ( 1 tr) for the complete set and is 0.24 dex for 
the Puis et al. set. This is an expected scatter and it im- 
plies that we do not find a systematic discrepancy between 
observations and our predictions for the O star mass-loss 
rates. 

Contrary to earlier comparisons between observations 
and theory where systematic discrepancies have been re- 
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ported (see Lamers & Leitherer 1993, Puis et al. 1996), 
here we find that there is agreement between our predic- 
tions and the mass-loss rates derived from observations for 
the O-type stars. The essential difference between previ- 
ous studies and the present one is that in our treatment of 
the theory of line driven winds, we consistently take into 
account effects of "multiple-scattering" in the transfer of 
momentum from the radiation field to the wind. We find 
systematic agreement between observed and theoretical 
mass-loss rates for a large sample of O stars. This result 
implies that physical effects that were not incorporated in 
our models, such as magnetic fields and stellar rotation, 
is not expected to influence the mass-loss rates of O stars 
significantly. 

6.2. Modified Wind momentum comparison for Range 1 
(27 500 < T cB < 50 000 K) 

Instead of comparing just the mass-loss rates it is useful 
to compare (modified) wind momenta derived from ob- 
servations and theory. In earlier studies, e.g. Lamers & 
Leitherer (1993), and Puis et al. (1996), wind momenta 
have been plotted versus the wind efficiency number r\. 
Comparisons between observed and theoretical wind mo- 
menta as a function of r\ could yield important informa- 
tion about the origin of the systematic discrepancy be- 
tween theory and observations. However, since these two 
quantities (wind momentum and wind efficiency number) 
both contain the mass-loss rate, they are not independent. 
Therefore, no such comparison is made here. Instead, the 
wind momenta are plotted versus the stellar luminosity, 
to compare the observational and theoretical WLR. 

We divide the T e s range into two parts. First, we ex- 
amine the wind momenta for stars where T e g > 27 500 
K, later on we will also compare the cooler stars. Fig- 
ure ^| shows the modified wind momentum as a function 
of stellar luminosity for the sample of stars with known ob- 
servational mass-loss rates. The upper panel shows these 
modified wind momentum values for the theoretical mass- 
loss rates and a linear best fit through these theoretical 
data (dotted line). Note that the "theoretical" WLR only 
contains the theoretical M, the included values for Voo and 
i?* were taken from observations. The theoretical WLR is: 
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Fig. 9. Upper panel: The theoretical modified wind mo- 
mentum expressed in Mq /yr km/s i?0 ' 5 for the stars in 
the first T eff range (27 500 < T eB < 50 000 K). The dotted 
line indicates the best linear fit. Lower panel: The obser- 
vational modified wind momentum for these stars. The 
dotted line indicates the same theoretical linear fit, as in 
the upper panel. 



jjthcory = _ 12 J2 ( ± q_ 26 ) 

+ 1.826 (± 0.044) log(L/L Q ) 

for T cff > 27 500 K (15) 

Since the slope of the WLR of Eq. ( |l5| ) has a slope of 
x = 1.826, the derived theoretical value for a c s (Eq. |l0|) 
that follows, is: 

a cS = - = 0.548 (16) 



This corresponds well to predicted values of the force mul- 
tiplier parameter (a ~ 0.66 and 5 ~ 0.10, see e.g. Paul- 
drach et al. 1994). 

The lower panel of Fig. ^| shows that both the WLR for 
the Puis et al. (1996) data and that for the other meth- 
ods/authors, follow the same relationship, both in agree- 
ment with the theoretical WLR. The dotted line is again 
the theoretical best linear fit. We conclude that for the 
range of the O stars, there is good agreement between 
theoretical wind momenta and those determined from ob- 
servations. The scatter between theoretical and observa- 
tional modified wind momenta is only 0.06 (1 a). 
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The good agreement between the observational and 
theoretical wind momenta adds support to the possibility 
to derive distances to luminous, hot stars in extragalactic 
stellar systems using the WLR. In practice the technique 
may be hampered by e.g. the fact that O stars are mostly 
seen in stellar clusters and cannot be spatially resolved 
in distant stellar systems. This is one of the reasons why 
the visually brighter B-type and especially the A-type su- 
pergiants located in the field are expected to be better 
candidates in actually using the WLR as a distance indi- 
cator (see Kudritzki et al. 1999). 

Comparison between the theoretical and observational 
WLR for the winds of B and A type supergiants is thus es- 
sential to investigate whether the slope of the WLR is the 
same for different spectral ranges. This is not expected, 
since the winds of different spectral types are driven by 
lines of different ions (see Vink et al. 1999; Puis et al. 
2000). 

6.3. Modified Wind momentum comparison for Range 2 
(12 500 < T cS < 22 500 K) 



Figure 10 shows the modified wind momentum as a func- 
tion of luminosity for both theory and observations for 
the stars in the second range (12 500 < T cff < 22 500 K). 
A best fit through the theoretically derived WLR is indi- 
cated with a dotted line in both panels. The theoretical 
WLR for this T c g range is: 



jjtheory = _ u 2g ( ± q 23) 

+ 1.914 (± 0.043) log(L/L ) 

for 12 500 < T cff < 22 500 K (17) 

Since the slope of the WLR for this range is slightly higher 
than that for the O star range, the predicted value for a e g 
is somewhat lower (see Sect. 4.3), namely: 



ttoff 



1 

X 



0.522 



(18) 



The lower panel of Fig. |0| indicates the observed mod- 
ified wind momenta (the dotted line contains the theoret- 
ical mass-loss rates). For this second T e s range (12 500 
< Tcff < 22 500 K) the plot in the lower panel reveals a 
large scatter in the observed data. 

Comparison of these observations with our predictions 
shows that within the subset of radio mass-loss rates there 
does not appear to be a systematic discrepancy. Also, 
those Ha profiles which are fully in emission (the filled 
symbols in the lower panel of Fig. [H]), i.e. the profiles 
that within the Ha method most likely provide the most 
reliable mass-loss rates, do not show a systematic differ- 
ence with the radio rates. The picture becomes different 
for stars that show Ha to be P Cygni shaped (grey symbols 
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Fig. 10. Upper panel: The theoretical modified wind mo- 
mentum expressed in M Q /yr km/s Rq ' 5 for the second 
range (12 500 < T eS < 22 500 K). The dotted line indicates 
the best linear fit for this range. The solid line indicates 
the theoretical WLR for stars in the range 27 500 < T e ff < 
50 000 K. Lower panel: The observational modified wind 
momentum for stars in this T c ff range (12 500 < T e ff < 
22 500 K). The different sources of the observations are 
indicated in the upper left corner. The dotted line again 
indicates the theoretical linear fit for the stars in the sec- 
ond range (12 500 < T cff < 22 500 K). 



in lower panel of Fig. |l0| ) or fully in absorption (open sym- 
bols). Although the measurements of Scuderi (1994,2000) 
remain reasonably consistent, those by Kudritzki et al. 
(1999) are discrepant in that at log L/Lq ~ 5.8 these val- 
ues start to diverge from the other observed rates, such 
that below log L/L Q ~ 5.6 a systematic difference of 
about a factor of 30 results between different sets of ob- 
served mass-loss rates. 
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An investigation of the origin of these systematic dif- 
ferences in observed B star wind momenta is beyond the 
scope of this paper. We will address this issue in a sep- 
arate study (Lamers et al., in preparation). Here we just 
note that the large scatter in the observed Ha data implies 
that there is either a dichotomy in the wind-momenta of 
B-stars (as suggested by Kudritzki et al. 1999) or that 
there exist systematic errors in the mass-loss determina- 
tions from Ha for B stars. 

The systematic discrepancies for the observed B star 
wind momenta imply that we cannot currently compare 
our predictions with observed data in the most meaning- 
ful way, since the data are not consistent and thus a fair 
comparison with our predictions cannot be conclusive. In 
addition, it may be meaningful to further investigate the 
validity of assumptions in our method of predicting the 
mass-loss rates of B-type stars (see e.g. Owocki & Puis 
1999). Still, we note that the most reliable rates (from ra- 
dio and pure Ha emission profiles) appear to be consistent 
with our predictions. 

The upper panel of Fig. [l(] reveals that most of the 
models in the second T cff range (12 500 < T off < 22 500 
K) lie above the theoretical fit for the models from the first 
T e g range (T c g > 27 500 K). This is due to the increase 
in the mass-loss rate at the bi-stability jump of a factor 
of five. The models with 12 500 < T eS < 22 500 K are, 
however, only slightly above the fit for the O star models 
(^eff > 27 500 K), as at the bi-stability jump the terminal 
velocity Voo drops by a factor of two. 

7. Discussion 

We have shown that our predictions of mass loss for O 
stars, using Monte Carlo simulations of energy loss during 
photon transport in non-LTE unified wind models, yields 
good agreement with the observed values. This demon- 
strates that an adequate treatment of "multiple scatter- 
ing" in radiation-driven wind models resolves the discrep- 
ancy between observations and theory that had been re- 
ported earlier. The agreement between observed and the- 
oretical wind momenta of O stars adds support to the 
method of deriving distances to distant stellar systems 
using the WLR. 

The comparison between the predicted and observed 
values of the modified wind momentum n for the B stars 
is not conclusive. A good comparison between the obser- 
vations and our predictions for the B star regime needs to 
await an explanation of the discrepancies in the observed 
B star mass-loss rates. This issue will be addressed in a 
separate study. 

Our models predict a jump in mass loss of about a 
factor of five around spectral type Bl. An important point 
that supports this prediction is the following. Vink et al. 
(1999) have calculated the mass-loss rate and Woo for winds 
at both sides of the bi-stability jump in a self-consistent 
way for models with log(L/i ) = 5.0 and M = 20 M Q . 



These self-consistent calculations showed a jump in M of 
a factor of five and a simultaneous drop in Voo/v csc of 
a factor of two. This drop in Woo/v CS c has been observed 
(Lamers et al. 1995). This gives support to our prediction 
that the mass-loss rate at spectral type Bl increases by 
the predicted amount. 

Since there is good agreement between observed mass- 
loss rates by different methods and the new theoretical 
mass-loss rates for the O-type stars, whereas there is in- 
consistency between the observed mass-loss rates from dif- 
ferent authors for the B-type stars, this may point to the 
presence of systematic errors in mass-loss determinations 
from observations for B stars. 

Because our predictions for the O stars agree with 
observations and our models also predict the bi-stability 
jump around spectral type Bl, we believe that our theo- 
retical mass-loss predictions are reliable and suggest they 
be used in new evolutionary calculations of massive stars. 

8. Summary &: Conclusions 

1. We have calculated a grid of wind models and mass- 
loss rates for a wide range of stellar parameters, cor- 
responding to masses between 15 and 120 Mq. 

2. We have derived two fitting formulae for the mass- 
loss rates in two ranges in T e g at either side of the 
bi-stability jump around 25 000 K. A mass-loss recipe 
was derived that connects the two fitting formulae at 
the bi-stability jump. 

3. There is good agreement between our mass-loss pre- 
dictions that take multiple scattering into account, and 
the observations for the O stars. There is no systematic 
difference between predicted and observed mass-loss 
rates. 

4. A comparison between observed and predicted wind 
momenta of O-type stars also shows there is good 
agreement. This adds support to the use of the WLR 
as a way to derive distances to luminous O stars in 
distant stellar systems. 

5. For the observed mass-loss rates of B stars there is an 
inconsistency between rates derived by different au- 
thors and/or methods. One group of M determina- 
tions of B stars does follow the theoretical relationship, 
while another group does not. This lack of agreement 
between the observed mass-loss rates of B stars may 
point to systematic errors in the observed values. 

6. Since our new theoretical mass-loss formalism is suc- 
cessful in explaining the observed mass-loss rates for 
O-type stars, as well as in predicting the location (in 
T e ff) and size (in Doo/fcsc) of the observed bi-stability 
jump, we believe that our predictions are reliable and 
suggest that our recipe be used in new evolutionary 
calculations for massive stars. A computer routine to 
calculate mass loss is available either upon request or 
at the following url: www.astro.uu.nl/~jvink/. 
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